AD-A142  544 
UNCLASSIFIED 


-  ■  -  11  k1 

THEORY  OF  THE  ELECTRONIC  STRUCTURE  AND  OPTICAL  1/1 

PROPERTIES  OF  ORGANIC  SOLIDS:  EXCITON  EFFECTS(U) 

MICHIGAN  TECHNOLOGICAL  UNIV  HOUGHTON  A  B  KUNZ 
30  JUN  84  N00014-81 -K-0620  F/G  7/4 


NL 


MICROCOPY  RISOLUMON  ItSt  CHARI 
NATIONAL  Hum  All  oi  SIANOMHS  “ 


\ 


A142  544 

*  ■  I  nr  #ri  •  i 


SCCuniT  V  CL  ASSIFIC  a  i  ION  Ol  Tins  r*M  'W**#*" 


REPORT  DOCUMENTATION  PAGE 


t.  REPORT  NUMBER 


■»  Title  t~>d  Submit) 


'heory  of  the  Electronic  Structure  and  Optical 
’roperties  of  Organic  Solids:  Exciton  Effects 


KKAO  INSTI  .'CHONS 
BEKokk  ruMi'i.n  in«~,  t-n 


*  TYPE  of  RtPOBT  4  PERIOD  COVERED 

Interim  Report  y 

1  July  1983  -  30  June  1984 


4.  PERFORMING  one.  REPORT  NUMBER 


author?*) 

1.  Barry  Kunz,  Presidential  Professor 
’hysics  Dept.,  Michigan  Technological  University 
loughton,  MI  49931 

PERFORMING  ORGANIZATION  name  and  aooress 

|  Michigan  Technological  University 
Houghton,  MI  49931 


»«.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Office  of  Naval  Research 
Arlington,  VA  22217 


4.  MONITORING  AGENCY  NAME  a  AOORESS(f/  dlllermnl  from  C  onlrulllttg  Ol/lrm) 


•  .  CONTRACT  OR  GRANT  MUM9E  Rf*) 


GNR-N-00014-81-K-0620 


10.  PROGRAM  ElCmCnI  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 


IZ  REPORT  OATE 

30  June  1984 


•  >  NUMBER  OF  PAGES 

17 


IS  SECURITY  CLASS,  (ol  Ihlt  upon) 


Unclassified 


is*,  decl  assifi  cation/  oowngr  aoing 

SCHEDULE 


14.  Distribution  statement  <oi  <M»  R.po ti> 


Approved  for  public  release:  distribution  unlimited. 


>7.  OlST  Rl  0  ;j  TIC  N  STATEMENT  (of  <hm  mbmtrmcl  mnirntmd  In  Block  20.  II  dlilmttnl  Itom  Rmpott) 


If.  KtY  »OROS  (Conf/noo  on  rommimm  mldm  It  nocmmtmry  mnd  Identity  by  block  numbor) 

Organic  Solids,  Methane,  excitons,  relaxation,  instability 


20.  ABSTRACT  ("Cofrtfow*  an  MftrM  mid*  H  nacaaaao'  Mrf  Identity  by  blurb  numbotf\ 


he  effect  of  exciton  formation  in  solid'1  "Hydro-Carbons  is  studied  and  structural 
instabilities  which  accompany  such  are  examined.  An  ablnitio  theory  of  such 
effects  is  presented. v-  ... 


I 


4f\  r*  O 

_ \j  ii  -»v 


11  3 


DO  ,',r *„  1<73  COITION  OF  •  NOV  4S  IS  OBSOL ETC 


security  classification  of  this  page  o»>» 


*in  concentractions  and 
ion  between  electrons  due  to 
the  exchange  of  virtual 
mb  corresponding  to  an 
an  between  the  electrons, 
for  the  superconductivity. 


liparticle  called  an  "exciton" 
ling  the  physics  and  chemistry 
jood  definition  to  some 
can  be  used  in  describing  the 
ir  the  superconductivity 


sincere  thanks  to  Professor 

ma. 


17. 

7)  191. 

■on”,  Reynolds,  O.C.  and 
c.,New  York)  1981. 

d). 

•J.»  Chandrasekhar,  M.  and 
27  (1983)  140. 

648. 


THEORY  of  the  electronic  structure  and  optical 

PROPERTIES  OF  ORGANIC  SOLIDS:  COLLECTIVE  EFFECTS* 


A.  Barry  Runs 

Department  of  Physics  and  Materials  Research  Laboratory 
University  of  Illinois  at  Urb ana-Champalgn 
Urbaaa,  Illinois  61801,  U.S.A. 

Abstract:  In  this  series  of  lectures  we  briefly  consider/. two 

complimentary  approaches  to  the  study  of  organic  solids:  The 

method  of  simulation  by  finite  clusters  of  molecules,  and  the 
methods  of  energy  band  theory.  In  both  cases,  the  initial 
starting  point  is  the  Hartree-Fock  method,  which,  as  expected, 
turns  out  to  be  inadequate  for  any  reasonable  level  of  quanti¬ 
tative  accuracy. n  Solids,  being  essentially  infinite  sized 
'Systems,  restrict  our  choice  of  correlation  methods  to  those 
i  which  are  size  consistent.  We  are  furthermore  interested  in 
(  properties  such  as  optical  excitation  and  need  to  be  able  to 
.obtain  the  finite  difference  between  extensive  total  energies. 
This  further  restricts  our  choices  .-^Methods  based  upon  ordinary 
Raylelgh-Schrodinger  Perturbation-Theory  are  chosen  and  exten¬ 
sive  results  for  solid  CH^  are  used  as  an  illustration.^— 

This  research  has  been  funded  in  part  by  the  U.  S.  Navy  Office 
•f  Naval  Research,  0NR-N-0014-81-K-0620,  in  cooperation  with  the 
Department  of  Physics,  Michigan  Technological  University,  and  by 
the  National  Science  Foundation,  DMR-80-20250  in  cooperation 
with  the  Materials  Research  Laboratory  of  the  University  of 
Illinois. 


I*  Introduction. 

Theoretical  studies  on  the  electronic  structure  of  three 
^■•nsional  solids  have  largely  excluded  the  organic  or  molec- 
■lnr  solids.  The  vast  majority  of  existing  calculations  have 
performed  for  the  solid  rare  g8ses.  More  complicated 
■oleeular  solids,  such  as  those  with  two  or  more  atoms  per  molecu- 
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lar  wit,  or  more  than  one  molecular  unit  per  unit  cell,  have 
been  largely  Ignored,  The  principal  exception  to  this  tendency 
baa  been  for  solid  It  should  be  further  noted,  that  the 

Interest  In  solid  H2  stems  largely  from  Interest  In  Its  possible 
transformation  to  a  monatomic  metal,  exhibiting  high  temperature 
BCS  type  superconductivity.  In  addition,  a  recent  resurgance 
of  Interest  in  the  solid  rare  gases  has  been  generated  by 
speculatlons  that  "low"  pressure  (~0.3  M  bar)  metal  phases  of 
solid  Xe  might  exist.  In  addition,  some  theoretical  studies  on 
solid  N,4  and  (^Cr  exist.  In  addition,  quite  a  few  studies  on 
properties  of  polymer  systems  exist.”  There  are  probably  many 
reasons  for  the  neglect  of  this  Interesting  and  technologically 
Interesting  class  of  solids.  Several  of  the  reasons  are  likely 
related  to  the  complicated  and  at  times  Ill-defined  crystal 
Structure  of  such  systems  and  the  associated  difficulties  In 
constructing  adequate  theoretical  models.  A  second  and  perhaps 
■ore  serious  problem  relates  to  the  question  of  which  approach 
one  might  use  to  determine  the  electrical  structure.  As  an 
example,  the  spectrum  of  solid  CH^  has  been  determined  over  an 
energy  range  of  8  to  about  35  eV.  The  fundamental  spectral 
region  of  from  threshold  (>8.5  eV)  to  about  14  eV  shows  marked 
Similarity  In  both  solid  and  gas  phase.  It  is  generally 
conceded  that  the  gas  phase  spectra  in  this  energy  region  Is 
dominated  by  transitions  from  the  bonding  to  antibonding  bound 
State  orbitals  or  to  Rydberg  series  like  transitions.  It  seems 
reasonable  to  expect  that  the  crystalline  spectrum  is  likely  to 
be  similarly  dominated  by  transitions  to  bound  rather  than  free 
final  states.  That  is  we  do  not  expect  the  contributions  from 
energy  band  theory  to  play  a  major  role  In  the  low  lying 
excitations  of  solid  CH^.  On  the  other  hand  the  spectral  region 
above  14  or  so  eV  may  well  be  dominated  by  band  to  band 
transitions  and  this  may  account  for  the  apparent  differences 
between  the  high  lying  spectrum  of  gas  phase  CH4  and  solid 
CHi.  Similar  considerations  apply  to  many  other  molecular 
solids. 


The  previous  theoretical  study  on  solid  methane  lends  cred¬ 
ibility  to  this  argument,  as  the  calculation  of  Plela, 
Pletronero,  and  Resta  finds  a  band  gap  In  excess  of  27.2  eV  for 
solid  methane.  It  is  not  likely  that  this  result  Is  quanti¬ 
tatively  accurate  as  these  authors  used  a  very  abbreviated  basis 
set  In  their  calculation  and  found  the  conduction  results  to  be 
highly  sensitive  to  the  virtual  basis  set.  A  further  source  of 
error  In  this  early  study  Is  the  use  of  the  Hartree-Fock  approx¬ 
imation  uncorrected  for  any  correlation  corrections.  Similar 
studies  by  Mlckish  and  Kunz  on  the  somewhat  similar  solid  rare 
gases  have  found  that  the  Hartree-Fock  method  consistently  over¬ 
estimates  the  band  gap  of  these  systems  by  about  4  or  5  eV.  In 
addition  all  band  methods  are  inaccurate.  In  that,  all  neglect 


A.B.KWZ 


COLLECTIVE  EFFECTS  IN  SOLIDS 


ts 


molt  cell,  hanre 
Do  this  tendency 
noted,  that  the 
:  In  Its  possible 
high  temperature 
■ecent  resurganee 
so  generated  by 
Metal  phases  of 
ideal  studies  on 
e  few  studies  on 
re  probably  many 
technologically 
neons  are  likely 
■defined  crystal 
difficulties  In 
eond  and  perhaps 
C  which  approach 
me t  tire .  As  an 
teralned  over  an 
■nntal  spectral 
off  shows  marked 
rt  Is  generally 
nargy  region  la 
sdbondlng  bound 
Urns. 1  It  seema 
Tai  Is  likely  to 
rather  than  free 
strlbutlona  fron 
the  low  lying 
spectral  region 
'  bend  to  band 
■wnt  differences 
CH4  and  solid 
other  molecular 


■■a  lends  cred— 
loo  of  Plela, 
of  27.2  eV  for 
■wit  le  quantl- 
■brwvlated  basis 
e  results  to  be 
rther  source  of 
'  no-Fock  approx— 
Ions.  Similar 
:  liar  solid  rare 
alstently  over- 
yl  or  3  mV.  In 
.  it,  all  neglect 


t 


f 


f 


i 


f 

i 

i 

» 

■ 


i 

? 

i 

i 

i 

i 

; 

l 

i 


i 

t 

* 

t 

* 


i 

» 


the  foraetlon  of  local  excited  states  called  excltons. 

He  believe  new  approaches  are  needed  If  one  Is  to  truly 
Interpret  the  electronic  structure  of  such  systems  as  solid 
Methane.  Recent  theoretical  results  of  Runs  and  Flynn  have 
deaonstrated  that  It  Is  possible  to  Include  the  effect  of 
electron-hole  Interaction  and  exclton  formation  without  violat¬ 
ing  Bloch's  theorem  in  calculations  of  the  optical  properties  of 
such  divergent  solids  as  L1F  and  Mg  or  Ca.  This  is  accomplished 
by  aeans  of  s  degenerate  perturbative  calculation  using  the  &- 
conserved  one  body  valence  to  conduction  band  excitations  as  a 
basis.  This  model  retains  the  periodic  symmetry  of  the  lattice 
avoiding  complications  Introduced  by  the  use  of  finite  cluster 
models  to  describe  the  local  excitations.  These  finite  cluster 
models,  nonetheless,  are  useful  and  accurate  approximates  as  we 
shall  see.  The  formulations)  of  the  problem  In  thl9  way  by  Kunz 
and  Flynn  causes  one  to  wi9h  to  begin  with  Hartree-Fock  descrip¬ 
tions  of  the  solid  since  a  well  defined  Many-Body  wavefunctlon 
Is  needed.  The  Hartree-Fock  model  neglects  all  correlations  and 
the  limited  basis  set  used  to  describe  excltonlc  effects  does 
not  describe  properly  the  relaxation  or  polarization  properties 
of  the  system.  In  these  lectures,  the  correlation  effects  are 
Incorporated  by  means  of  a  simple  Many  Body  Perturbation  Theory 
Model  (MBPT).  The  necessary  theoretical  methods  are  described 
In  Section  II.  The  numerical  calculations  are  described  in 
Section  III,  and  conclusions  are  given  In  the  final  section. 

II.  Theoretical  Development 

The  Initial  step  In  this  development  Is  the  choice  of  the 
Hartree-Fock  method.  This  choice  Is  largely  determined  by  the 
need  to  perform  extensive  correlation  calculations  In  addition 
to  Che  Initial  Hartree-Fock  study.  To  facilitate  development, 
wo  employ  variants  on  the  familiar  Llnear-Comblnatlon-of-Atomlc- 
Orbitals  method  (LCAO).  In  the  case  of  cluster  calculations, 
these  AO's  are  first  rotated  Into  molecular  orbitals  (MO's) 
■panning  the  entire  cluster,  and  in  the  case  of  the  band  calcu¬ 
lations,  the  AO's  are  rotated  Into  MO's  spanning  the  crystallo¬ 
graphic  unit  cell.  ThiB  rotation  Into  MO’s  Is  advantageous  be¬ 
cause  for  unit  cells  of  ever  Increasing  size  or  complexity,  an 
adequate  description  In  terms  of  AO's  yields  rather  large 
secular  determinants.  The  LCMO  scheme  reduces  substantially  the 
•1m  of  the  secular  determinant.  This  method  was  first  Intro¬ 
duced  by  Plela  _et_  al.  for  studies  on  solid  methane.  In  such  a 
simple  case  the  basis  set  for  the  occupied  orbitals  Is  reduced 
from  9  to  S  orbitals.  Furthermore,  the  MO's  may  contain  polari¬ 
zation  functions  In  them  and  therefore  yield  far  greater 
accuracy  than  a  much  larger  set  of  AO's. 
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The  essential  features  of  this  approach  is  this.  Let  each 
unit  cell  be  divided  into  molecules  (real  ones  or  not),  and  we 
devise  a  basis  set  to  represent  the  MO 's  of  these  molecules. 
The  primitive  basis  set  used  are  spherical-harmonic  Gaussian 
Type  Orbitals  (GTO's)  centered  about  different  origins,  and  have 
the  form 

-  exp  [-Z^Cr-fti)2]  Y*®  (6,*).  (1) 


The  Hi  are  the  origins  about  which  these  functions  are  centered 
end  sued  not  be  an  actual  nuclear  site,  the  Y,m  are  the  usual 
spherical  harmonics.  The  orbital  exponent  2^  is  chosen  by 
energy  minimization.  The  MO's  in  turn  are  just  linear 
combinations  of  these  GTO's, 

♦j  (r-Hj,)  -  £  a^  (2) 

The  H^'s  are  the  locations  of  the  molecules  in  the  system.  Prom 
the  MO's,  one  forms  Bloch  orbitals  which  span  the  entire  system: 

tj(H,r)  -  IT1'2  l  e1***0  (3) 

The  MO's  (Eq.  (2))  or  the  Bloch  functions  (Eq.  (3))  form  the 

basis  by  which  we  solve  the  HF  problem  or  its  extensions. 

The  first  point  is  that  the  Hartree-Fock  equation  need  be 
solved  self-conslstently.  For  a  finite  molecular  cluster,  this 
is  achieved  by  conventional  iterative  means.  However,  the 

infinite  periodic  system  imposes  special  difficulties.  These 

are  simply  that  the  occupied  canonical  Bloch  orbitals  are 

infinite  in  number  and  therefore  enumerating  the  contribution  of 

sach  orbital  to  the  Pock  operator  Imposes  a  strain  on  ones  com¬ 
puter  budget.  Two  options  are  available.  The  first  is  to  use  a 
finite  mesh  in  k-  space  and  use  some  form  of  quadrature  to 

construct  the  Fock  operator.  The  second  is  to  rotate  into  a 

basis  set  of  local  orbitals.10,11’12  In  the  early  stages,  both 
methods  were  tried  with  negligible  differences  in  numerical 

result  between  them.  However,  at  the  current  stage  of  our  code 
development,  the  local  orbitals  method  enjoys  a  substantial 

speed  advantage. 

t 

The  intent  of  the  present  study  is  to  obtain  spectroscopic 
information  and  hence  we  need  examine  the  meaning  of  the  energy 
bands.  The  occupied  bands  are  the  negative  of  the  ionization 
energy  for  that  band  for  the  state  of  wave  vector  £  .  The 

virtual  bands  are  similar  representations  for  the  electron 

affinities.  In  this  event  one  is  assuming  the  use  of  the 

Hartree-Fock  eigenvalue  and  also  of  Koopmsns'  theorem  as  is 
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usually  done.  The  essential  physics  here  refers  to  Ionization 
properties,  not  to  excitation  properties  of  the  n-electron 

system. 

In  order  to  improve  upon  the  Hartree-Fock  results  one  must 
Include  correlation  corrections.  In  doing  this.  Initially  the 
author  will  maintain  the  same  physical  definition  for  the  energy 
bands  as  In  the  Koopmans'  case.  That  Is,  the  bands  now  become 
quasl-partlcle  bands  In  which  Che  energy  of  an  occupied  level  Is 
the  negative  of  Che  energy  needed  to  create  It,  and  the  energy 
of  the  virtual  states  are  the  negative  of  the  energy  recovered 
In  creating  It.  This  Is  in  keeping  with  the  earlier  usage  of 
the  electronic  polaron  model  and  Its  extensions  as  discussed  by 
Pantelldes  et  al. 13,14,15 
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It  Is  now  necessary  to  discuss  correlation  corrections. 
The  first  problem  is  that  of  .lire  consistency  (The  total  energy 
Is  an  extensive  quantity).  1 '  In  fact  the  total  energy  of  an 
infinite  solid  Is  infinite  and  only  the  energy /molecule  Is 
finite.  Unfortunately,  the  energy  change  upon  ionization  Is 
also  finite  and  the  energy  change/molecule  vanishes.'  That  Is 
the  energy  difference  Is  still  finite.  Similar  considerations 
apply  to  excitations  of  the  n-electron  system.  A  simple  clas¬ 
sical  way  to  view  this  is  to  realize  that  the  size  of  the  wave 
created  by  hurling  a  brick  Into  a  pond  Is  largely  Independent  of 
the  size  of  the  pond.  Therefore  we  must  establish  a  size  con¬ 
sistent  framework  for  the  system  total  energies  In  such  a  way 
that  formally  we  can  obtain  differences  In  extensive  quantities, 
cancelling  the  Infinities  before  we  compute  finite  differences. 
Alternately,  we  may  reduce  the  size  of  the  systems  so  that  total 
energy  determinations  are  possible. 


Let  us  work  In  a  local  representation  here.  This  Is  appro¬ 
priate  since  many  molecular  solids  are  filled  shell  systems. 
For  notational  simpllcitly,  designate  the  Uannler  function 
W.jjCr)  as  the  1th  Wannler  function  about  site  Form  a  com¬ 

plete  set  of  Wannler  orbitals  describing  the  ground  state  of  the 
neutral,  N-electron  solid  in  the  Hartree-Fock  limit.  Me  will 
use  them  to  generate  the  Ion  states  as  well.  For  a  system  of  N- 
electrona  the  Hamiltonian  Is 


f  -iiyj-f  f  J-I--  +  i  f  f  1  _al_  (4) 
1-1  “  1  1-1  1-1  |r-dj|  *1-1  j-1  |?jj  | 


The  electronic  has  mass  m,  and  Is  chars 


Is  the  atomic 


*•  - 

number  of  the  nucleus  at  site  I.  The  l*-  electron  has  coordin¬ 
ate  r^  and  the  I*1*  nucleus  has  coordinate  ftj.  In  terms  of 

Vannler  Functions,  In  the  single  determinant  limit,  the  energy 
of  the  system  Is 
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eh  -JP  <  -ill-  £  72  -  X  !-«>  '  (5> 

+  i  ^  f  ^ii-jjIf^l-iiHj^  -  <  -iiwjjlfYjlwjj«u>  1 

The  symbol  (N)  on  the  summation  implies  sums  over  all  states  in 
the  occupied  N  electron  space.  To  keep  the  physics  of  the  en¬ 
ergy  bands  discussed  earlier  we  need  to  look  at  the  N-l  and  N+l 
electron  system  next. 

Let  the  ground  state  of  the  N-electron  Hartree-Fock  system 
be  designated  as  |N>  and  let  a+  »,  a_j  create  or  destroy  a 
Wannler  function  at  site  I  with  other  quantum  numbers  p.  We 
adopt  the  conventlve  that  quantum  numbers  i,  j,  k  etc.,  refer  to 
occupied  orbitals,  a,  b,  c  to  virtual  orbitals  and  o,  p,  q  to 
either/both.  A  Slater  determined  of  the  N-l  body  system  is 

|R-1;  JB  >  =  ajB|N  >  ’  (6) 

This  will  by  symmetry  adapted  later.  The  energy  expectation 
value  of  this  state  is  simply 

Here  F(N)  is  simply  the  N-electron  ground  state  Hartree-Fock 
operator.  Similarly  one  may  obtain  the  off  diagonal  matrix 
elements  between  two  states  |N-1,  1A>  and  |N-1,  JB>.  These  are: 

“  <wJB|F(N)|w1A>.  (8) 

One  aiay  project  on  the  state  |N-1.  lb>  to  form  a  proper  transla¬ 
tional  invariant  Bloch  function,  ♦j  ^  (1c) : 

<*)  -  J  ^R*i*’*B  |W-1;  JB>  (9) 

In  terms  of  eqs  (5)-(9)  one  may  construct  a  band  structure  in 
terms  of  Wannler-functions  and  Slater  determlnante  for  the 
occupied  orbitals.  These  are  yet  uncorrelated.  One  may  treat 
the  N+l  body  states  similarly.  Furthermore,  recognizing  that  E,, 
la  eq  (11)  is  Infinite  and  also  Irrelevant,  since  energy  changes 
ote  needed,  we  proceed  to  define  EQ  as  0,  and  thus  simplify  our 
computation.  1 

A  framework  is  needed  in  order  to  simply  correlate  this 
problem  since  the  simple  single  Slater  determinants  |n-l,  JB> 
are  highly  degenerate,  and  within  a  band,  the  p,(()  are  nearly 
degenerate.  Consider  the  problem  in  a  general  framework 
Initially.  H  Is  a  Hamiltonian, 
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H  -  H  +  V.  (10) 

o 

lie  assume  that  the  eigenstates  of  HQ  are  known  as 

Vi  “  Vi*  (I1) 

The  projector  onto  a  given  eigenstate  of  HQ  say  4^  is  P^  and  la 
given  as 

pi  “  iV^i1* 

Furthermore  a  projector  onto  the  first  n  states  say  is  P  and 


n  n 

p  -  l  lv<*il  •  l  pr 

i-1  11  i-1 


Assume  we  order  our  eigenfunction  of  Hq  so  that  the  states  of 
Interest  lie  In  the  range  0  to  n.  Furthermore  no  other  states 
are  degenerate  with  these  states.  Now  let  us  solve  the  desired 
equation: 

fKr  -  E*  -  (H0+V)<|i.  (I3) 

Let  us  chose  a  Wj  for  1  <  j  <  n,  then: 

<l-?)(H  -w,)*  -  (l-P)(E-v.-V)<p  (14) 

o  J  J 

One  may  commute  (1-P)  with  H0-Wj  ar,d  proceed  to  see 

t  -  P*  +  (H#-m  )-1  (1-P)  (E-Wj-V)*.  (15) 

Furthermore; 

-  «1f1. 

so  that 


P*  -  I  Vk  “ 

k-I 


Therefore 


[l  -  (l#-w  )-1  (l-PKEitj-P)]*  ” 


If  one  defines 


T  -  (1  -  (H  -w  J-1  ( l-P)(E-w.-V))_l , 


jtcpr.-- 
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than 


♦  -  t*. 


Furthermore  one  can  show  that 
n 


<E"V  "  l  "k  Vtv* 

x  i  k-1  K 


where 


09) 

(20) 


«*> 

Eqs.  (20)  and  (21)  define  a  perfectly  good  algebraic  eigenvalue 
equation  for  the  system  energies.  To  proceed  further,  one 
expands  the  Inverse  appearing  In  T.  That  Is, 

VT  -  N  +  V  (l-P)(E-w  -V)  +  ...  .  (22) 

o  J  J 

-  ?. 

or 


lk 


m 

*  l 

a-fC-fl 


— =—  V, 
Wj-wa  la 


ak* 


(23) 


The  structure  of  the  eigenvalue  problem  defined  by  Eqs  (20), 
(21)  and  (23)  Is  now  clear.  The  matrix  elements  to  lowest 
approximation  are  similar  to  those  of  second  order  R.S.P.T.  and 
this  Is  clearly  a  size  consistent  approach.  If  all  the  eigen¬ 
vectors  In  the  first  n  are  degenerate,  one  recovers  normal  de¬ 
generate  perturbation  theory.  Consider  our  problem,  where  we 
use  Wannler  functions,  this  framework  makes  our  approach 
clear.  First  correlate  the  single  Slater  determinant  of  Wannler 
functions,  then  proceed  with  Bloch  symmetry  projection  to  remove 
the  degeneracy.  The  N-body  wavefunctlon  has  proper  Bloch 
symmetry  for  closed  band  systems.  By  using  a  proper  choice  of  A 
in  the  Adams-Gllbert  local  orbital  formulation  called  Aw  one  may 
obtain  Wannler  orbitals.11  The  actual  choice  of  Aw  Is  not 
Important,  only  that  such  exist.  Then 

IF  +  PAWP]  wtI  -  etIwj|  (24) 

The  first  order  Fock-Dlrac  density  matrix  Is  p.  From  this  one 
constructs  a  zero  order  Hamiltonian.  For  a  system  of  M- 
electrona,  HQ  Is  defined  as 


Ho  ‘  J*  +  ’iW 


l 

V 

j 


(25) 
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gad  then  the  perturbation,  V  becomes 


?  5  H  -  H 


,aAbB  1 2 


(N)  |V“T 

EOO-E  +  I  [  e  .  (27) 

"  1I>JJ  aA>bB  HjJaA  aB 

Here  the  summation  (N)  Indicates  all  Wannler  orbitals  In  the  N- 
electron  state.  The  matrix  element  Is  simply 

^5?  ■  ^ilV^1  WaAWbB>-<WlIWJJl7^|wbBWaA>  (28) 

The  N-body  orbitals  will  be  used  to  describe  both  the  N-l  and 
H+l  body  states.  Brlllouin’s  theorem  Is  not  valid  for  such 
states.  Consider  first  the  N-I  body  problem.  Let  Wannler 
orbital  w,B  be  deleted  from  the  N-body  ground  state.  Then  to 
second  order  one  finds  that  .  7 

Of+D  r  |F(N-1;1B).J 

E(H-l.lB)  -  E  +  I  I  - r~=c’J - 1 

® ,  1  JI  aA  jl  aA 


(H-l) 

+  l 


kK>JJ  aA>bC  ekK+ejJ_eaA~ebC 

In  Eq  (19)  the  V  Is  still  as  defined  in  Eq  (18)  and  F(N-1;1B)  Is 
obtained  by  deleting  terms  referring  to  orbital  w^g  ;  rom  F(N). 
Therefore 

P(I»-1;1B)^  -  <wjI|F(N-l,lB)|waA>.  (30) 

One  proceeds  In  like  fashion  for  the  NH  electron  case,  adding 
W~  to  the  N  electron  state.  .  , 

b  b  OH-l)  1 F(N+1 ,cB)?j| 2 
E(IH-l;cB)  -  EN  +  +  l  [  - - 


JJ  aA 


JJ  aA 


OH-D  r 

♦  l  l 


JJ>W  aA>dD  EjJ4€kK~€aA"eaD 

la  Eq  (21)  V  remains  as  in  Eq  (18)  and  F(N+l;cB)  is  obtained  by 
adding  terms  referring  to  orbital  wcB  to  the  N-electron  Fock 
operator. 

One  may  obtain  the  physically  Interesting  energy 
differences  from  these  expressions.  The  Ionization  potentials 
•re  defined  by  E(N)  -  E(N-1;1B).  This  difference  called 
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here  4..  Is  given  as 
**  /» 


niB«  <Nf1}  v 

N_1  jJ  eA  EJJ  ~  caA 


1^0,2  i_iBaAi2 

♦y  i  iSi il _ r  i  IZimL 

JJ*1B  aA>cC  ClB'f€jJ_£aA~ecC  J^fgK  aA  ej J^kK^iB^aA 

Likewise  the  electron  affinity  terms  are  obtained  by  letting 
A  -  E(N+1;CB)-E(N).  Then 

rRr.  (N+l)  | F(N+1 ,  CB)“;r  (N+l ) 

‘a-C*  l  I  ■-  JJ...  .  i  i 

jJ  aA  jJ  aA  11  aA>dD*cB 


AcB  "  Vl 


(N+l) 

-  I  l 


,cBaA  1 2 


*11  +  *cB  “  £aA  “  edD  1I>JJ  aA*cB  E1I  +  EjJ  "  EcB  '  EaA 
It  is  these  formulas  we  will  use  In  this  study. 

One  final  piece  Is  needed  to  complete  this  theory.  This  Is 
to  include  the  actual  effect  of  electron-hole  interaction  upon 
excitation.  An  accurate  method  of  doing  this  for  both  tightly 
bound  or  loosely  bound  excitations  has  been  recently  given  by 
Kunz  and  Flynn. 

The  essential  point  Is  to  use  the  Hartree-Fock  bands  as  a 
basis  set  after  incorporation  of  correlation  corrections  into 
the  hand  energies.  The  Fock  ground  state  |N>  is  then  used  to 
describe  schematically  the  process.  Let  av(£)  annihilate  a 
valence  electron  of  wavevector  tc  and  let  a£(tc)  create  a  conduc¬ 
tion  electron  of  wavevector  &.  Consider  the  state  then: 

|H,t>  -  «+<*)  ay(lc>  |N>  (34) 

It  Is  only  states  like  this  which  can  be  reached  from  the  ground 
state  via  optical  processes.  Furthermore  all  such  ground 
states  |N,ti>  correspond  to  the  same  total  crystal  wavevector; 
that  of  the  ground  state.  The  most  general  excited  state  that 
one  nay  access  Is  then  |N,  E>,  where 


!*.«>  -  j. 


In  this  sum,  the  ground  state  (N>  Is  excluded  because  It  differs 
In  parity  from  the  excited  state.  By  finding  the  a£  and  <N. 
I|H|N,  E>,  one  may  determine  the  spectrum  of  the  solid  Including 
electron-hole  Interaction.  This  Is  achieved  by  means  of  a  Cl 
calculation  among  the  states  |N,  fc>.  The  formation  of  such  ex- 
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clton  states  is  not  an  extensive  property  and  size  consistency 
is  not  a  problem  as  demonstrated  by  Kunz  and  Flynn.  Exact 
Implementation  of  such  an  Infinite  Cl  is  of  course  Impossible 
and  we  use  a  finite  number  of  states,  some  270  configurations. 
A  second  approximation  Is  made  as  well.  This  Is  to  truncate  the 
coulomb  Interaction  at  the  boundary  of  a  unit  cell.  This  is  not 
unreasonable  for  tightly  bound  excited  systems  as  in  the  case  of 
CH^  particularly  since  the  large  lattice  constant  (11.14  au)  en¬ 
closes  a  substantial  volume  In  a  unit  cell.  The  dominant  conse¬ 
quence  of  this  is  to  allow  the  formation  of  only  a  single  bound 
exclton,  not  an  entire  Rydberg  series  below  the  bands.  However 
when  the  coefficients  a£  are  used  to  evaluate  the  optical  re¬ 
sponse  one  finds  substantial  adjustment  over  the  Hartree-Fock 
results.  These  changes  are  due  to  the  redistribution  of 
oscillator  strength  to  the  bottom  of  the  conduction  band  due  to 
the  Inclusion  of  electron-hole  interaction. 

The  alternate  approach  employed  is  to  use  a  finite  molecu¬ 
lar  cluster  simulation.  This  is  also  done  using  the  method  of 
local  orbitals.  In  this  case  one  partitions  the  system  into  the 
cluster  and  the  environment.  The  environment  imposes  itself  on 
the  cluster  by  means  of  a  bounding  potential.  The  methods  of 
doing  this  are  well  represented  in  the  literature,  and  a  general 
approach  is  given  by  Kunz  and  Kleinla  which  need  not  be  repeated 
here.  Correlation  is  imposed  using  the  technique  of  this  sec¬ 
tion  and  in  particular  eqs.  (20),  (21)  and  (23)  as  needed.  For 
non-degenerate  states  of  course,  these  reduce  to  ordinary  second 
order  RSPT.  Most  cases  considered  here  are  not  degenerate  in 
the  cluster  limit,  but  for  those  cases  for  which  degeneracy  is  a 
problem,  we  have  found  the  full  approach  to  be  very  powerful.1 


III.  Results  for  Solid  CH^ 

A  Gaussian  basis  set  was  first  developed  For  the  CH^  mole¬ 
cule  in  free  space  and  then  reoptimized  for  the  crystal  to  allow 
accurate  description  of  the  energy  bands,  occupied  and  virtual. 
It  was  found  easy  to  obtain  accurate  valence  bands,  but  that  the 
conduction  bands  were  quite  sensitive  to  the  choice  of  outer  or¬ 
bital.  The  variational  principal  applies  to  the  solution  of  the 
one  particle  states  in  a  LCMO  formalism,  and  the  selection  of 
the  basis  is  quite  easy.  In  practice,  the  conduction  bands  are 
found  to  be  stable  against  small  changes  in  basis  set.  The  va¬ 
lence  structure  here  agrees  well  with  .that  obtained  after 
corrections  to  formalism  by  Plela  j!t_  al .  *20,21  The  conduction 
bands  are  in  very  poor  agreement  however.  This  is  due  to  the 
far  coo  restrictive  basis  set  employed  in  the  Plela  et  al.  cal¬ 
culation  of  the  virtual  bands.  In  performing  this  calculation, 
•oae  Idealizations  are  needed.  A  lattice  constant  of  11.14  au. 


*4 


A.B.KUNZ 


In  agreement  with  Plela  _et_  al.  Is  used  and  the  C  sub  lattice  is 
fixed  as  a  fee  one  as  per  experiment.  The  four  H's  form  In 
tetrahedra  about  the  C  in  a  unit  cell.  In  the  real  world,  the 
tetrahedra  do  not  allign  from  one  cell  to  another  but  have 
orientational  disorder.  We,  as  did  Piela  et  al.  fix  the  H's  in 
an  ordered  fee  lattice  as  well..  The  current  calculation  uses 
the  same  geometry  as  does  Plela  et  al.  The  C-H  distance  is  ob¬ 
tained  computationally  from  Beck^^- and  for  a  lattice  constant  of 
11.14  au,  the  equilibrium  constant,  is  essentially  the  same  C-H 
distance  as  in  the  free  molecule. 

Although  the  Hartree-Fock  band  results  overestimate  any 
reasonable  band  gap,  they  do  reduce  the  Piela  gap  by  about  13.6 
eV  however,  and  one  need  add  correlation  along  the  lines 
suggested  in  Section  II.  In  performing  the  correlation  correc¬ 
tion  computations,  the  author  deviates  from  the  ideals  expressed 
In  the  preceding  section  to  the  extent  that  Instead  of  solving 
for  a  set  of  rather  complicated,  orthogonal  Wannier  functions  as 
implied  by  the  derivations,  one  approximates  these  by  a  set  of 

local  orbitals.  In  obtaining  these  the  unit  on  which  localiza¬ 
tion  occurs  is  the  CH^  molecule  is  used,  and  also  the  appropri¬ 
ate  multicenter  localization.  i  These  orbitals  are  quite  local, 
the  valence-valence  overlaps  being  0.03  or  less  here.  First 
order  overlap  corrections  are  made  in  the  inter  molecular  terns 
for  further  precision.  Due  to  the  procedure  adopted,  all  orders 
of  overlap  in  the  large  intra  molecular  overlaps  are  Included 
exactly.  The  inclusion  of  these  corrections  is  essential  if  one 
wishes  to  achieve  quantitative  accuracy.  In  evaluating  the 

perturbation  sums,  d  orbitals  on  the  C  atom  and  p  orbitals  on 
the  H  atoms  were  added  to  the  band  structure  basis  set.  The 
effect  of  the  several  contributions  to  Eqs  (32)  and  (33)  are 
given  In  Table  1.  There  we  call  the  second  term  on  the  right 
hand  sides  of  Eqs  (32)  and  (33)  the  relaxation  and  the  sum  of 

the  second  and  third  terras,  which  come  from  two  electron  virtual 

excitations,  clearly  represent  correlation  terms. 

The  energy  bands  for  CH^  including  correlation  are  shown  in 
Plgure  1.  The  density  of  electron  states  is  also  seen  in  these 
figures.  As  is  clear  from  these  figures,  the  band  gap  is  in¬ 
direct  and  from  to  Xj,,..  The  direct  gap  is  at  the  X  point 

and  is  Xj,y  to  Xj,£.  The  correlated  indirect  gap  is  found  to  be 
13.0  eV,  and  the  correlated  direct  gap  is  found  to  be  13.3  eV. 

Finally,  one  computes  the  position  of  the  exclton  levels  in 
CH^.  This  is  accomplished  using  the  method  given  in  Section  II 
which  has  been  more  fully  described  in  Ref.  9.  In  this  calcu¬ 
lation  the  coulomb  interactive  is  treated  as  a  one  molecule 
interaction.  The  effective  electron-hole  interaction  is  here 
computed  to  be  S.4  eV.  This  is  the  value  of  the  VQ  discussed  in 
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Table  l 

Contributions  to  the  ionization  potential  and  electron 
affinity  of  solid  CH^  are  shown  as  a  function  of  lattice 
parameter.  Results  are  given  for  the  correlation  correction  and 
the  relaxation  correction.  Results  are  in  eV. 


■ 

11.14  au 

10.50  au  |  10.00  au 

valence  correlation 

0.1  eV 

0.2  eV  |  0.4  eV 

valence  relaxation 

1.2  eV 

1.2  eV  |  1.2  eV 

conduction  correlation 

-.7  eV 

-.8  eV  |  -.9  eV 

conduction  relaxation 

~0.0  eV 

~0.0  eV  j  ~0.0  eV 

net  gap  change 

-2.0  eV 

-2.2  eV  |  -2.5  eV 

Figure  1  The  correlated  band  structure  of  solid  CH^  and 

density  of  states  in  shown  for  lattice  parameter  = 
11.14  au. 
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I ef.  9.  Using  this  value,  Che  exclton  Is  found  to  be  at  10.9  eV. 
The  optical  spectrum  of  Koch  and  Skihowski24  does  find  a 
epectral  peak  at  11.0  eV  and  this  may  well  be  our  exclton.  A 
■ore  quantitative  analysis  of  the  optical  response  Is  not 
possible  at  this  time  because  Che  highest  valence  and  lower 
conduction  bands  are  of  like  symmetry  and  the  techniques 
developed  in  Ref.  9  and  currently  available  do  not  permit  a  cal¬ 
culation  of  the  optical  response  for  the  case  in  which  the  band 
Co  hand  transitions  are  dipole  forbidden  only  the  positions. 
Therefore,  the  author  reluctantly  contents  himself  with  using 
only  the  k  conserved  Joint  density  of  states  in  comparison  with 
the  eeasured  reflectance  spectrum  shown  in  Figure  2. 24  As  Is 
clear  from  this  figure,  even  if  one  were  to  Include  the  exclton 
et  11.0  eV,  a  fair  degree  of  discrepancy  remains.  This  is 
largely  at  low  energy.  This  discrepancy  is  expected.  A  similar 
result  is  seen  in  the  free  CH^  molecule  and  is  due  to  the 
■oblllty  of  the  H  nuclei  and  their  large  zero  point  motion.  The 
excited  CH^  molecule  can  lower  its  energy  by  about  1.6  eV  by 
relaxing  from  ideal  TQ  geometry  to  Djj,  geometry.  Due  to  large 
xero  point  motion  it  may  be  possible  to  excite  from  the  ground 
state  TD  geometry  directly  into  the  relaxed,  distorted  D,h 
geometry  directly.  This  certainly  appears  to  be  the  case  in  the 
free  molecule  and  a  discussion  of  this  is  being  prepared  by  Beck 
and  Kunz.  5  If  one  assumes  the  same  type  of  Jahn-Teller  distor¬ 
tion  is  present  in  the  solid,  a  distorted  exclton  line  would 
then  appear  at  about  9.3  eV.  This  is  shown  as  a  dotted  line  in 
Figure  i.  Since  the  first  experimental  peak  in  solid  CH^  lies 
at  9.6  eV,  this  inclusion  greatly  enhances  the  comparison  of 
theory  and  experiment.  In  addition  the  low  energy  continuous 
apectrum  between  about  12  and  14  eV  would  be  enhanced  in 
strength  by  the  redistribution  of  oscillator  strength  due  to 
exclton  formation  as  was  seen  in  L1F. 

Large  scale  cluster  calculation  for  bulk  CH^  (13  molecules 
or  65  atoms)  and  for  the  CH^  surface  (9  molecules  or  45  atoms). 
Including  all  electrons  and  correlation  via  the  perturbative 
roots,  have  been  recently  performed  by  Beck.2"  These  calcu¬ 
lations  are  for  the  excitons  alone  and  tend  to  confirm  the 
energy  band  results  qualitatively  and  quantitatively.  The 

specific  details  of  the  perturbation  treatment  for  large  systems 
la  sell  described  in  the  literature. 27 
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IV.  Conclusion* 

The  essential  conclusions  are  few  and  simple.  These  are 
•ae  can  construct  a  satisfactory,  self-consistent  Hartree-Fock 
hand  structure  for  molecular  solids,  including  the  conduction 
heads,  if  one  carefully  optimizes  the  basis  set.  If  one  wishes 
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to  obtain  quantitative  comparisons  with  experiment,  the 
Inclusion  of  correlation  corrections  Is  essential.  Furthermore, 
in  describing  the  ion  states  In  terms  of  the  neutral  system 
orbitals  corrections  termed  relaxation  corrections  are  needed. 
It  Is  seen  here,  using  a  Wannler  basis,  how  such  arise  and  may 
be  included.  It  Is  also  seen  that  Inclusion  of  electron-hole 
Interaction  Is  needed  If  one  Is  to  quantitatively  study  the 
optical  spectrum.  In  addition,  due  to  the  light  mas  of  H  one 
need  also  be  prepared  to  Include  Jahn-Teller  distortion  If  one 
Is  to  be  fully  quantitative. 
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Figure  2  The  optical  joint  density  of  state  is  shown  for  solid 
Cfy  along  with  the  T0  geometry  exciton  position  and 
probable  D2n  geometry  exciton.  The  optical  reflectivity 
of  Ref.  24  is  also  shown 
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